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Summary. Lattice QCD predicts a phase transition between hadronic matter and 
a system of deconfined quarks and gluons (the Quark Gluon Plasma) at high en- 
ergy densities. Recent results from the Brookhaven Relativistic Heavy Ion Collider 
(RHIC) dedicated to the study of QCD at extreme densities will be discussed and 
compared to measurements obtained at the CERN Super Proton Synchrotron (SPS). 



1 Introduction 

Quantum Chromo Dynamics (QCD) provides, as part of the standard model, 
a very successful description of strong interaction processes involving large 
momentum transfer. However, from first principles several important aspects 
of QCD are still poorly understood. Examples are color confinement, chiral 
symmetry restoration and the structure of the vacuum. Better understanding 
of these concepts can be obtained if we are able to study quarks and gluons 
in a deconfined state, the so-called Quark Gluon Plasma (QGP). 

Such a deconfined state might be created in the laboratory in heavy-ion 
collisions at the highest energies. Theoretical guidance for this comes from 
Lattice QCD calculations. Lattice QCD predicts that at an energy density 
e « 1 GeV/im^, corresponding to a temperature of about 170 MeV, the system 
undergoes a phase transition from nuclear matter to a deconfined system of 
quarks and gluons. 

Figure^ shows the energy density divided by the fourth power of the tem- 
perature, T, versus the temperature from a lattice calculation Jj. This figure 
shows that in between 150-200 MeV the energy density increases rapidly which 
is indicative of a phase transition where at high temperature the quarks and 
gluons become the relevant degrees of freedom. The figure also indicates where 
according to our current understanding the different heavy-ion machines are 
located on this diagram. These calculations are done with zero baryon chemi- 
cal potential, hb, reflecting the conditions of the early universe. Small values 
of the chemical potential are obtained at RHIC collider energies whereas at 
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Fig. 1. a) Energy density divided by T* versus T at /ib = d'- b) Transition 
temperature as a function of The dotted line illustrates the rapid crossover 

while the solid line illustrates the first order phase transition. 



lower energies, e.g. AGS and SPS, the value of fiB is large. In Fig. ^ the 
relation between the transition temperature and the chemical potential from 
recent lattice calculations |21 is shown. The calculation indicates that the tran- 
sition temperature decreases with increasing fj,B and furthermore that at low 
He the transition from the hadronic phase to the QGP is a rapid crossover 
(dotted line) while at large a first order transition should take place (full 
line). 

2 The Relativistic Heavy Ion Collider (RHIC) 

Heavy-ion physics entered a new era with the advent of the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory. RHIC is a versatile 
collider providing collisions with different ion species (ranging from protons 
to gold) at a wide range of center of mass energies y/s^. In the four years of 
operation collisions were provided for Au-|-Au at 19.7, 130 and 200 GeV, p-l-p 
at 200 GeV and d-l-Au at 200 GeV. Note that the top center of mass energy 
for p-l-p is 500 GeV at RHIC. For details on the machine and the detectors 
(BRAHMS, PHENIX, PHOBOS, STAR), see |S||lj. 

3 Event Characterization 

Heavy ions are extended objects and the system created in a head-on colli- 
sion is different from that in a peripheral collision. Therefore, collisions are 
categorized by their centrality. Theoretically the centrality is characterized by 
the impact parameter b which, however, is not a direct observable. Experi- 
mentally, the collision centrality can be inferred from the measured particle 
multiplicities if one assumes that this multiplicity is a monotonic function of 
b (see Fig.|2K)- Another way to determine the event centrality is to measure 
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the energy carried by the spectator nucleons (which do not participate in the 
reaction) with Zero Degree Calorimetry (ZDC), shown in Fig. ^p. A large 
(small) signal in the ZDCs thus indicates a peripheral (central) collision. 




Fig. 2. a) Multiplicity distribution measured in STAR The different colors de- 
note the different fractions of the cross section, b) ZDC signal versus multiplicity, 
measured by PHENIX 



Two other measures of the centrality which are often used are the number 
of wounded nucleons and the equivalent number of binary collisions. These 
numbers are related to the impact parameter b using a realistic description of 
the nuclear geometry in a Glauber calculation, see Fig. 13 Phenomenologically 
it is found that soft particle production scales with the number of participating 
nucleons whereas hard processes scale with the number of binary collisions. 



4 Low-pt Observables 

Examples of global observables which provide important information about 
the created system are the particle multiplicity and the transverse energy. 
Figure 0] shows the transverse energy versus the collision centrality as mea- 
sured at ^/s^ = 130 GeV by the PHENIX collaboration 0. This measure- 
ment allows for an estimate of the energy density as proposed by Bjorken 

_ 1 1 dEr 
■kE? ctq dy 

were R is the nuclear radius and tq is the effective thermalization time (0.2-1.0 
fm/c). From the measured {dET /drj) = 503 ± 2 GeV it follows that e is about 

5 GeV/fm'^ at RHIC. This is much larger than the critical energy density of 
1 GeV/fm3 from Lattice QCD (see Fig.[TJ. 

Figure |S1 shows the charged particle multiplicity distributions versus the 
pseudorapidity 77 measured by PHOBOS at three different energies [H]. The 
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Fig. 3. Number of wounded nucleons Fig. 4. Transverse energy as a func- 
and binary collisions versus impact pa- tion of centrality as measured by 
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Fig. 5. Multiplicity versus pseudo-rapidity for 19.6, 130 and 200 GeV measured by 
PHOBOS ^. 



gross features of the particle multiplicity distributions are described by a 
similar behavior of the tails (limiting fragmentation) and a plateau at mid- 
rapidity consistent with a boost invariant region of Ay k, 1. Notice that in 
total about 5000 charged particles are produced in the most central Au+Au 
collisions at the top RHIC energy. 

4.1 Particle Yields 

The integrated yield of the various particle species provides information on 
the production mechanism and the subsequent inelastic collisions. A very suc- 
cessful description of the relative particle yields is given by the thermal model. 
In Fig.|n|the particle yield ratios measured at RHIC are plotted and compared 
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Fig. 6. Particle yield ratios at RHIC compared with a thermal model 



to values from a thermal model fit JI] . The results from the fit show that all 
particles ratios are consistent with a single temperature and single chemical 
potential in a thermal description. The temperature obtained in this way, 176 
MeV, is called the chemical freeze-out temperature. Figure [3l shows the rela- 
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Fig. 7. a) Particle ratios versus rapidity measured by BRAHMS }1(J( - b) The ratio 
K~ /K~^ versus p/p or, equivalently, fis- 



tive particle ratios of pions, kaons and protons and their anti-particles versus 
rapidity JD]. For the protons and kaons the ratio drops rapidly for y > 1. 
Figure shows the ratio of K~ / versus p/p for AGS to RHIC energies. 
The decreasing ratio of p/p as a function of rapidity can thus be understood 
from the changing baryon chemical potential at a constant chemical freeze-out 
temperature. 
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compared to the critical temperature from lattice QCD calculations. 



Figure |S1 shows the chemical freeze-out temperature obtained in the ther- 
mal model versus in the SIS to RHIC energie range. The chemical freeze- 
out temperature increases strongly from SIS to SPS energies above which it 
seems to saturate close to the phase boundary temperature from lattice calcu- 
lations (see also Fig.QJj). This observation is not inconsistent with the scenario 
that the matter produced at SPS and RHIC energies was first thermalized in 
the deconfined quark-gluon plasma phase and subsequently expanded through 
the phase boundary into a thermal gas of hadrons. For a detailed overview of 
particle production and the thermal model see |13j . 

4.2 Spectra 

The particle spectra provide much more information than the integrated par- 
ticle yields alone. The particle yield as a function of transverse momentum 
reveal the dynamics of the collision, characterized by the temperature and 
transverse flow velocity of the system at kinetic freeze-out. Kinetic freeze-out 
corresponds to the final stage of the collision when the system becomes so 
dilute that all interactions between the particles cease to exist so that the 
momentum distributions do not change anymore. Figure and b show the 
transverse momentum distributions at -^s = 17 GeV from NA49 ^j. The 
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Fig. 9. NA49 (SPS) and RHIC low-pt spectra fT^I fT?^ fl^ 171 f7^ [Hj . 



lines are a fit to the particle spectra with a hydrodynamically inspired model 
(blast wave) . The fit describes all the particle spectra rather well which shows 
that these spectra can be characterized by the two parameters of the model: a 
single kinetic freeze-out temperature and a common transverse flow velocity. 
FigureEl; shows the combined pion, kaon and proton pf-spectra from the four 
RHIC experiments. Also at these energies it follows from a common fit to all 
the spectra that the system seems to freeze-out with a similar temperature 
and transverse flow velocity as observed at SPS energies. 



4.3 Azimuthal Correlations with the Reaction Plane 

The nuclear overlap region, shown in gray in Fig. llOh .. has in non-central col- 
lisions an almond like shape with its longer axis perpendicular to the reaction 
plane (the plane defined by the beam axis Z, and the impact parameter b). 
This particular shape leads to a pressure gradient which is different in and out 
of the reaction plane which, in turn leads to azimuthally asymmetric particle 
emission. The asymmetry can be described by: 

where (j)' is the azimuthal angle with respect to the reaction plane and the 
coefhcient of the second harmonic, V2, is called elliptic flow. The magnitude 
of V2 and its pt dependence allows for the extraction of the kinetic freeze-out 
temperature and the transverse flow velocity as function of emission angle. 

Figure llOb shows the integrated value of V2 versus beam energy. The 
magnitude of V2 increases smoothly from AGS to the top RHIC energy. At 
the highest RHIC energies for the first time in heavy-ion collisions the value 
reaches qualitative agreement with prediction from hydrodynamic model cal- 
culations 21 . 

Figure 1111 shows the measurement of V2 versus pt for pions and protons 
plus antiprotons. Due to transverse flow the pt dependence of V2 depends on 
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Fig. 12. V2 (pt ) for identified particles 
compared to a blast wave fit \2'A \24\ i25[ 



the particle mass as is evident from Fig. Also shown in this figure are 
hydrodynamical model calculations using two different equations of state |22| 
corresponding a hadron gas and a QGP. It is seen that the QGP EOS shows the 
best agreement with the data. In Fig.[l3RHIC data on W2(]3t) ISHIEIEHIEEI 
for various particles are compared to a hydrodynamical inspired blast wave 
fit. The agreement of the data with this fit shows that the f2(pt) for all 
particles can be described in terms of a single temperature and a (/)-dependent 
transverse fiow velocity. Furthermore, the magnitude and pt dependence of 
the elliptic flow for the various particles suggest strong partonic interactions 
in an early stage of the collision and, perhaps, early thermalization of the 
system. 
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4.4 Hanbury-Brown Twiss Interferometry 

Two-particle intensity interferometry, or the Hanbury-Brown Twiss |27) (HBT) 
effect, is a technique used to measure the size of an object emitting bosons. In 
heavy-ion collisions pion HBT has been used extensively to probe the space- 
time structure of the produced system. 




For a given pair of pions, we can define their momentum difference, q, and 
their momentum average, k (see Fig. 113b'). With these two vectors and the 
beam direction we can define the coordinate system used in HBT [211 : The 
longitudinal direction (i?;), which is along the beam direction z. The outward 
direction {Ro), in the z,k plane and _L z. Finally the side- ward direction 
(i?s), _L z and _L k. At low-pt, for a boost-invariant source, the side- ward 
radius will correspond to the actual physical transverse (RMS) size of the 
source at kinetic freeze-out. At larger-pt the source size reflects the region of 
homogeneity, due to the transverse flow of the system. The outward radius 
contains a mixture of the spatial and time extent of the source. Figure 113b 
shows the energy dependence of the tt~ HBT radii [2H1- The evolution of 
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the HBT radii shows a smooth dependence versus center of mass energy. The 
observed ratio Rq/Rs is part of the so-called "RHIC HBT puzzle" . The value is 
very close to unity, which naively implies an almost instantaneous emission of 
particles. The models, which are successful in describing the measured spectra 
and elliptic flow, predict larger values for this ratio. However in a blast wave 
description the dependence of the HBT radii versus pt is consistent with the 
large transverse flow deduced from the identified particle spectra and V2 (23 • 

5 High-j?t observables 

In heavy-ion collisions at RHIC, jets with transverse energies above 40 GeV 
are produced in abundance, providing a detailed probe of the created system. 
However the abundant soft particle production in heavy-ion collisions tends to 
obscure the characteristic jet structures. At sufficient high-pf the contribution 
from the tails of the soft particle production becomes negligible and jets can 
be identified by their leading particles. It was proposed that a leading particle 
traversing a dense system would lose energy by induced gluon radiation (so 
called jet-quenching |3H). The amount of energy loss is in this picture directly 
related to the parton density (mainly gluons at RHIC) of the created system. 
Currently there are three observables sensitive to this energy loss as discussed 
in the next two subsections. 

5.1 Single Inclusive Particle Yields 

As mentioned above, the single inclusive particle yield at sufficiently high-pt is 
dominated by the leading particles from jets. Figure IT^ shows the ttq spectra 
as measured in p-|-p at ^/s — 200 GeV. In the same figure also two NLO QCD 
calculations are shown. The ratio of the data to the theory shows that in p+p 
the TTo spectrum is well described. In Fig. 114b the charged hadron spectra 
measured in Au-|-Au at ^s„„ = 200 GeV and the p-|-p reference spectra at 
the same energy are shown. One of the observables suggested for measuring 
energy loss is the so called nuclear modification factor defined by 

i?AA(Pt) = d^o-AA/d^/dpf 

(A^binary)d2(Tpp/dydpt ' 

where a-p-p / dydpt is the inclusive cross section measured in p4-p collisions 
(see Fig.lTlk.e') and (A^binary) (see Fig. EI accounts for the geometrical scaling 
from p-l-p to nuclear collisions. In the case that a Au-|-Au collision is an inco- 
herent superposition of p-l-p collisions this ratio i?AA would be unity. Energy 
loss and shadowing would reduce this ratio below unity while anti-shadowing 
and the Cronin effect would lead to a value above unity. Figure [TSb .c shows 
this ratio for charged particles and tto's in central Au-|-Au collisions at mid- 
rapidity. The ratio is well below one and at high-pt the suppression is a factor 
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Fig. 14. a-d) ttq spectra in p+p yj'^ compared to NLO QCD calculations. e)Charged 
hadrons in p+p and Au+Au ]3cl^ . 



of 5. At intermediate pt the charged particles and ttq are both suppressed; 
however the magnitude differs by a factor of two. In Fig. [TM l Raa is plotted 
at more forward rapidities showing that the suppression also persists there. 

To discriminate between energy loss and shadowing, d+Au collisions were 
measured. If the suppression is due to shadowing it should also be observed 
in the d-Au system. Figure 115b shows the d+Au spectra versus centrality 
and Fig. llSb.c the nuclear modification factor for charged particles and ttq, 
respectively. It is clear that in d+Au interactions no suppression is observed. 
In fact, to the contrary, a small enhancement is seen consistent with the 
Cronin effect. From this observation it follows that the observed suppression 
in Au+Au collisions is due to final state interactions. The magnitude of the 
observed suppression at the top RHIC energy indicates, in the jet quenching 
picture, densities which are a factor 30 higher than in nuclear matter. 

5.2 Azimuthal Correlations 

In heavy-ion collisions, azimuthal correlations between particles can be used to 
study the effect of jet quenching in greater detail. The azimuthal correlations 
of two high-p4 particles from jets are expected to show a narrow near-side 
correlation and a broader away-side correlation. However, in the case of strong 
jet quenching the away-side jet would suffer significant energy loss and would 
be suppressed. Recently, CERES measured such a correlation function at the 
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Fig. 16. Azimuthal correlations measured at the SPS by CERES f3^ . 



top SPS energy. In Fig. 116b the nearside correlation (at Acf) = 0) shows a 
narrow peak consistent with the correlation observed in jets. The away-side 
correlation peak is observed in more peripheral collisions but disappears for 
more central collisions, see Fig. 116b. FigurelTBb shows that the width (cr) of the 
near-side correlation peak stays constant as a function of centrality, but that 
the away-side peak broadens for more central collisions. The total integrated 
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yield is the same in the near and away-side peak (Fig. llGti'). Therefore, the 
disappearance of the away-side peak at the top SPS energy is interpreted as 
being due to initial state broadening |37j. 




Fig. 17. a,b) Back to back correlations jmilf and c) elliptic flow parameter ^241 at 
intermediate pt . 



The azimuthal correlations of high-pt particles (trigger particle 4 < pf < 
6 GeV/c, associated particle 2 GeV/c < pt < pl"^) measured in p-l-p collisions 
at RHIC are shown as the histogram in Fig. 117b. The near-side and away-side 
peaks are clearly visible. The correlation function observed in central Au+Au 
collisions (stars in Fig. 117b) shows a similar near side peak while the away-side 
peak has disappeared. 

To investigate if this is due to initial state effects, the same analysis was 
done for d-|-Au collisions. In Fig. 117b the near and away-side peaks are shown 
for minimum bias and central d-f-Au collisions compared to p-l-p. The away- 
side correlation in d-f-Au is clearly observed even for the most central col- 
lisions. Comparing the away-side correlation in p-l-p, d-|-Au and Au-|-Au, 
Fig. [T7b . shows that the suppression only occurs in Au-|-Au collisions and 
therefore is a final state effect as expected from jet quenching. 

The energy loss depends on the distance traversed through the dense 
medium by the partons. In a non-central collision the distance will depend 
on the azimuthal angle with respect to the reaction plane [221 (see \ow-pt sec- 
tion). Because the hard scattering producing the di-jet has no correlation with 
the reaction plane, an observed asymmetry in the high-pt particle emission 
will be due to final state interactions (such as the jet quenching mechanism). 
In Fig. [T7b the observed elliptic flow signal as a function of pt is shown for 
charged particles, kaons and lambdas. It is clear from this figure that the 
observed asymmetry is very large up to the highest pt measured. Like the 
nuclear suppression factor Raa, the elliptic flow at intermediate pt depends 
on the particle species. This could be due to an interplay between the soft 
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hydrodynamical behavior and the jet quenching, which would cause a mass 
dependence |4(Jj . However, more recently this has been interpreted as a possi- 
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Fig. 18. a) Charged pion, charged kaon and proton plus antiproton V2 measured by 
PHENIX 126L h) and lambda plus antilambda V2 measured by STAR '24^- Shown 
is V2 versus pt where both are normalized by the number of constituent quarks. 



ble sign of particle production at intermediate pt by parton coalescence |41| . 
In that case it is not the mass of the particle which is responsible for the 
splitting but rather the number of constituent quarks (two for mesons and 
three for baryons). Figure [THk.b shows the elliptic flow versus pt, where both 
are normalized by the number of constituent quarks [201 1^- Plotted like this, 
at intermediate pt, the V2 of the different species should reduce to an ap- 
proximately universal curve. The measurement of the pions, kaons, protons 
and lambdas V2 indeed seem to follow this scaling. A definitive test will be 
the measurement of elliptic flow of the 0-meson because in the coalescence 
interpretation it should have an elliptic flow similar to the pions while in the 
hydrodynamical interpretation is would have an elliptic flow value similar to 
the proton. 



6 Conclusions 

The first four years of RHIC operation have provided a wealth of interesting 
data. We have seen that: 

• Particle yields indicate a chemical freeze-out of the system near the phase 
boundary; 

• Identified particle spectra are consistent with boosted thermal distribu- 
tions and identified particle elliptic fiow shows remarkable agreement with 
ideal hydrodynamical calculations based on a QGP equation of state; 
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• The particle yield at high-p^ is suppressed compared to proton-proton 
reference data. The fact that this suppression does not occur in d+Au 
collisions shows that it is a final state effect, consistent with parton energy 
loss in dense matter (jet quenching); 

• The suppression at intermediate pt shows a particle dependence which 
could be explained by particle production, at intermediate pt, by parton 
coalescence; 

• The elliptic flow at intermediate pt is large and also shows a particle de- 
pendence. Like above, this is consistent with energy loss in dense matter 
and particle production via parton coalescence; 

• In the most central events the high-pt back to back correlations are con- 
sistent with zero. Such disappearance of the away-side jets is expected in 
the case of very strong energy loss in a dense medium. 

All these observations, taken together, are consistent with the creation of 
a very dense and strongly interacting system in heavy-ion collisions at RHIC 
energies. While all these observations are consistent with the creation of a 
QGP, more detailed knowledge of QCD at high densities and temperatures is 
required. This poses a formidable challenge for theory but will be crucial for 
the detailed interpretation of the present and future data taken at RHIC and 
LHC. 
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